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We report work on the matching of the next-to-leading logarithmic approximation (NLLA) onto the fixed next-to- 
next-to-lcading order (NNLO) calculation for event shape variables in electron-positron annihilation. The correction 
of the combined NLLA+NNLO computation in the three-jet region, relevant for precision phenomenology, is small 
compared with pure NNLO or NLLA+NLO. 

1. INTRODUCTION 

Event shape distributions in e + e~ annihilation processes are very popular hadronic observables, mainly due to the 
fact that they are well suited both for experimental measurement and for theoretical description because many of them 
are infrared and collinear safe. The main idea behind event shapes variables is to parameterize the energy-momentum 
flow of an event, such that one can smoothly describe its shape passing from pencil-like two-jet configurations, which 
are a limiting case in event shapes, up to multijet final states. Since the deviation from two-jet configurations 
is proportional to the strong coupling constant a s , the comparison of experimental measurements and theoretical 
prediction permits to determine a s . 

At LEP a set of six different event shape observables were measured in great detail: thrust T (which is substituted 
here by r = 1 — T), heavy jet mass p, wide and total jet broadening Bw and Bt, C-parameter and two-to-three-jet 
transition parameter in the Durham algorithm V3. 

Until very recently, the theoretical state-of-the-art description of event shape distributions was based on the 
matching of the NLLA Q onto the NLO 0,0] calculation. Using the newly available results of the NNLO 1 corrections 
for the standard set of event shapes 0] introduced above, we computed the matching of the resummed NLLA onto 
the fixed order NNLO. 

2. FIXED ORDER AND RESUMMED CALCULATIONS 

At NNLO the integrated cross section 

R{y,Q,v) = 

0"had 

has the following fixed-order expansion: 

R (y, Q, m) = 1 + « s (m) a (y) + (m) b (v, x n) + (m) c (y> • 

where a s = a s /(2n) and = (J,/Q. Approaching the two-jet region the infrared logarithms in the coefficient 
functions becomes large, spoiling the convergence of the perturbation expansion. The main contribution in this case 
comes from the highest power of the logarithms which have to be resummed to all orders. For suitable observables 
resummation leads to exponentiation. At NLLA the resummed expression is 




1 Recently an inconsistency in the treatment of large-angle soft radiation was discovered Q. It is about to be corrected and it should 
result in numerically minor changes to the NNLO coefficients in the kinematical region of phenomenological studies here. The corrections 
turn out to be significant only in the deep two-jet region, e.g. (1 — T) <C 0.05. 
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Table I: Powers of the logarithms present at different orders in perturbation theory. The color highlights the different orders in 
resummation: LL (red) and NLL (blue). The terms in green are contained in the LL and NLL contributions and exponentiate 
trivially with them. 



where the function g\ (a s L) contains all leading-logarithms (LL), g 2 (ot s L) all ncxt-to-leading-logarithms (NLL) and 
/i = Q is used. Terms beyond NLL have been consistently omitted. The resummation functions g\{a s L) and g 2 (a s L) 
can be expanded as power series in a s L 

L gi {a s L) = G 12 L 2 a s + G 23 L 3 o? s + G 34 L 4 a 3 s + . . . (LL) , 
g 2 (a s L) = G ll La s +G 22 L 2 a 2 s +G 33 L 3 a 3 s + ... (NLL). (1) 

Table U shows the logarithmic terms present up to the third order in perturbation theory. At the fixed order level 
the LL are term of the form a"L n+1 , the NLL those which goes like a™L n , and so on. Notice that this can be read 
off the expansion ([I]) of the exponentiated resummation functions. 

Closed analytic forms for functions qi (a s L), q 2 {a s L) are available for r and p 0], By/ and Bt [1, 0], C [h| and 



11 1, and are collected in the appendix of [l4j. Recently also g 3 (a s L) and 54 (a s L) were computed for r using 



effective field theory methods [3] 



3. MATCHING OF FIXED ORDER AND RESUMMED CALCULATIONS 

To obtain a reliable description of the event shape distributions over a wide range in y, it is mandatory to combine 
fixed order and resummed predictions. The two predictions have to be matched in a way that avoids the double 
counting of terms present in both. At NLLA the the expression which has to be matched with fixed NNLO is given 
by 

R(y)= (l + Cia s +<J 2 a;:+C3c^)e Lff1 ^^^^ D (y) , 

= C(a s )S(y) +D(y) , 

where S(a s ) is the exponentiated part containing the resummed logarithms, C (a s ) is a constant and D (cv s ) is a 
remainder functions which tends to zero as y — * 0. 

A number of different matching procedures have been proposed in the literature, see for example [l[ for a review. 
In the so-called i?-matching scheme, the two expression for R(y) are matched. In this case all the coefficients (C 2 , 
G3, G21, G32 and G31) appearing in @ have to be extracted numerically from the distributions at fixed order. The 
increasing number of logarithms present in the fixed order coefficient functions of the NNLO distributions causes 
large errors on these coefficients. For this reason we computed the matching in the so-called In i?-matching since 
in this particular scheme, all matching coefficients can be extracted analytically from the resummed calculation. The 
In i?-matching at NLO is described in detail in 0]. In the In i?-matching scheme, the NLLA+NNLO expression is 



In (R(y,a s )) = L g x (a s L) + g 2 (a s L) + a s (A{y) - G U L - G 12 L 2 

1 

2" 

+ a| (c (y) -A(y)B (y) + ^A 3 (y) - G 33 L 3 - G M L* 



a 2 s (B (y) - -A 2 (y) - G 22 L 2 - G 23 L 3 



(3) 



2 
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The matching coefficients appearing in this expression can be obtained from ([I]) and are listed in 14j . To ensure the 
vanishing of the matched expression at the kinematical boundary y max a further shift of the logarithm is made [if. 
The renormalisation scale dependence of ([3]) is given by making the following replacements: 

a s — > a s (fi) , 

B(y) -> B(y,n) = 20 o hix ll A(v)+B(y) , 

C(y) -> C(y,f,) = (2/3 ]nx fl ) 2 A(y) + 2]nx tl [2 /3 B (y) + 2 ft A (y)] + C (y) , 

g 2 (a s L) -> g 2 (a s L,fi 2 ) = g 2 (a s L) + — (a s L) 2 g\ (a s L) Inx^ , 

G22 -> G 22 (M) = G22 + 2/3 Gi2 ln Xfl , 

G33 G 33 (fi) = G33 + 4/3 G 23 In Xf, . 

In the above, g[ denotes the derivative of g\ with respect to its argument. The LO coefficient A and the LL 
resummation function gi, as well as the matching coefficients Ga + \ remain independent on [i. 



4. MATCHED DISTRIBUTIONS AND DISCUSSION 



For the resulting plots of the matched distributions we refer to [l4|. The most striking observation is that the 
difference between NLLA+NNLO and NNLO is largely restricted to the two-jet region, while NLLA+NLO and NLO 
differ in normalisation throughout the full kinematical range. This behavior may serve as a first indication for the 
numerical smallness of corrections beyond NNLO in the three-jet region. In the approach to the two-jet region, 
the NLLA+NLO and NLLA+NNLO predictions agree by construction, since the matching suppresses any fixed 
order terms. Although not so visible on these plots, the difference between NLLA+NNLO and NLLA+NLO is only 
moderate in the three-jet region. The renormalisation scale uncertainty in the three-jet region is reduced by 20-40% 
between NLLA+NLO and NLLA+NNLO. This effect is due to the smaller renormalization scale dependence of the 
NNLO contributions. It is also important to observe that the scale dependence remains the same and is larger in 
the two-jet region, because the resummed calculations at NLLA take into account only the one-loop running of the 
coupling constant. This has important consequences in the determination of a s and we will comment more on this 
in the next section. 

The description of the hadron-level data improves between parton-level NLLA+NLO and parton-level 
NLLA+NNLO, especially in the three-jet region. The behavior in the two-jet region is described better by the 
resummed predictions than by the fixed order NNLO, although the agreement is far from perfect. This discrepancy 
can in part be attributed to missing higher order logarithmic corrections and in part to non-perturbative corrections, 
which become large in the approach to the two-jet limit. 



5. CONCLUSIONS AND OUTLOOK 



After the extraction of a s using only the NNLO distributions and the experimental data of ALEPH [12||, a new 
extraction of a s using the new matched results was performed [ltl ] using JADE data. The improvement in the 
error coming from the inclusion of resummed calculation is not as drammatic as passing from NLO to NLLA+NLO 
calculations. As already anticipated, this is due to the fact that the NNLO coefficients compensate the two-loop 
renormalization scale variation, whereas the NLLA part only compensates the one-loop variation. A further im- 
provement is possible by including the NNLL corrections into the calculations. These corrections are known only for 
r, where higher order logarithmic corrections have been computed [l3| using soft-collincar effective theory (SCET). 
From these calculations one can extract the functions g% (a s L) and 174 (a s L). The next step towards the further 
improvement in the extraction of a s from event-shape distributions could be to compute them for all six observables 



mentioned here. As shown in I3J the subleading logarithmic corrections can also account for about half of the 



discrepancy between parton-level theoretical predictions and hadron-level experimental data. 
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Improvements can also come from non-perturbative corrections. A very recent non-perturbative study for r 



using a low-scale effective coupling 15j shows that non-perturbative 1/Q power corrections cause a shift in the 
distributions, which can account for an important part of the difference between parton-level distributions and 
hadron-level experimental data discussed in the previous section. 
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